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We consider slow-roll inflation for a single scalar field with an arbitrary potential and an arbitrary 
nonminimal coupling to the Gauss-Bonnet term. By introducing a combined hierarchy of Hubble 
and Gauss-Bonnet flow functions, we analytically derive the power spectra of scalar and tensor 
perturbations. The standard consistency relation between the tensor-to-scalar ratio and the spectral 
index of tensor perturbations is broken. We apply this formalism to a specific model with a monomial 
potential and an inverse monomial Gauss-Bonnet coupling and constrain it by the 7-year Wilkinson 
Microwave Anisotropy Probe data. The Gauss-Bonnet term with a positive (or negative) coupling 
may lead to a reduction (or enhancement) of the tensor-to-scalar ratio and hence may revive the 
quartic potential ruled out by recent cosmological data. 


PACS numbers: 98.80.Cq, 98.80.Jk, 04.62.+v 


I. INTRODUCTION 


Inflation in the early Universe has become the stan- 
dard model for the generation of cosmological perturba- 
tions in the Universe, the seeds for large-scale structure 
and temperature anisotropies of the cosmic microwave 
background. The simplest scenario of cosmological in- 
flation is based upon a single, minimally coupled scalar 
field with a flat potential. Quantum fluctuations of this 
inflaton field give rise to an almost scale-invariant power 
spectrum of isentropic perturbations (see Refs. [1, 2] for 
reviews). 


String theory is often regarded as the leading can- 
didate for unifying gravity with the other fundamental 
forces and for a quantum theory of gravity. It is known 
that the effective supergravity action from superstrings 
induces correction terms of higher order in the curvature, 
which may play a significant role in the early Universe. 
The simplest such correction is the Gauss-Bonnet (GB) 
term in the low-energy effective action of the heterotic 
string [3]. Such a term provides the possibility of avoiding 
the initial singularity of the Universe [4]. In the presence 
of an exponential potential for the modulus field, nonsin- 
gular cosmological solutions were found which begin in 
an asymptotically flat region, undergo superexponential 
inflation and end with a graceful exit to a phase with 
decreasing Hubble radius [5]. 

There are many works discussing accelerating cosmol- 
ogy with the GB correction in four and higher dimen- 
sions [6-9]. Recently it has been shown that the GB 
term might give rise to violent instabilities of tensor per- 
turbations [10]. A model in which inflation is driven by 
the GB term and a higher-order kinetic energy term was 
studied. When the GB term dominates the dynamics of 
the background, tensor perturbations exhibit violent neg- 
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ative instabilities around a de Sitter background on small 
scales, in spite of the fact that scale-invariant scalar per- 
turbations can be achieved [10]. Besides the kinetic and 
GB terms, a scalar potential arises naturally from super- 
symmetry breaking or other nonperturbative effects. 


In a previous work, we investigated inflationary so- 
lutions and resulting cosmological perturbations for the 
special case of power-law inflation when both the GB cor- 
rection and the scalar potential are present [11]. Power- 
law inflation happens when both the potential and the 
GB coupling take an exponential form. In this model 
instabilities of either scalar or tensor perturbations show 
up on small scales for GB-dominated inflation. The GB 
correction with a positive (or negative) coupling may lead 
to a reduction (or enhancement) of the tensor-to-scalar 
ratio in the potential-dominated case. This effect leads 
to tight constraints on the magnitude of the GB cor- 
rection from the Wilkinson Microwave Anisotropy Probe 
(WMAP) 5-year analysis [12]. 

Here we generalize our previous work to the more gen- 
eral case of slow-roll inflation with an arbitrary potential 
and an arbitrary coupling. Making use of a combined 
hierarchy (e€;, 6;) of Hubble and GB flow functions (as 
defined below) with |e;| < 1 and |é;| < 1, analogous 
to the standard slow-roll approximation, we derive the 
power spectra of scalar and tensor perturbations. In this 
scenario the spectral index of scalar perturbations con- 
tains not only the Hubble flow parameters but also the 
GB flow parameters. Moreover, the standard consistency 
relation of single-field slow-roll inflation is modified. In 
order to impose observational constraints on such models, 
we focus on a specific model with a monomial potential 
and an inverse monomial GB coupling. We analyze the 
influence of the GB term on the scalar spectral index ngr 
and the tensor-to-scalar ratio r. 

This paper is organized as follows. In Sec. II we define 
the Hubble and GB flow functions. Then by using the 
background equations of motion, we demonstrate that 
the slow-roll solution exists and is stable under asymp- 
totic conditions. In Sec. III we calculate the power spec- 


tra of scalar and tensor perturbations for the slow-roll 
inflation. In Sec. IV our approach is applied to a specific 
example. Section V is devoted to conclusions. 


II. SLOW-ROLL INFLATION WITH THE GB 
CORRECTION 


We consider the following action 
4 1 = 2 1 2 
S= | da/-g gift a av?) —V(¢) - zE) Rap (1) 


where ¢ is a scalar field with a potential V(¢), w = 
+1, R denotes the Ricci scalar, R43 = Ryvpo R” — 
4R uu R!” + R? is the GB term, and €(¢) is the GB cou- 
pling. We work in Planckian units, A = c = 87G = 1. 
In a spatially flat Friedmann-Robertson- Walker universe 
with scale factor a, the background equations read 


6H? = wọ? + 2V + 24ER”, (2) 
2H = —w¢? + 4€H? + 44H (2H = H?) i (3) 
i (è + 3H6) +Vg + 126 gH? (i + g) = 0, (4) 


where a dot represents the time derivative, (...),ẹ denotes 
a derivative with respect to ¢, and H = a/a denotes the 
expansion rate. Since the GB coupling is a function of ¢, 
one has € = £.4¢ and £ = E p° + Ego. 

Besides the slow-roll conditions ¢2 < V and |¢| « 
3H|¢|, well known for minimal-coupled single-field infla- 
tion, we impose two extra conditions, namely 4|€|H <1 
and || < |€|H. The background equations are approxi- 
mately given as 


1 
H? ~ -V 5 
Ly, (5) 
He -ug — 2¿R?, (6) 
; 1 E 
er -zg Me +12 gH”), (7) 


which allows us to obtain the number of e-folds 
$ 3wV 

N(¢) > 1 ae T 

bena 3V0 + 4E, V 


Following Ref. [13] we define a hierarchy of Hubble flow 
parameters, 


dọ. (8) 


H dln |e; | 
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21. (9) 
The expansion is accelerated as long as e€} < 1. In the 
slow-roll approximation they can be related to the usual 
slow-roll parameters. The new degrees of freedom intro- 
duced by the GB coupling function €(¢) suggest to define 


an additional hierarchy of flow parameters in the same 
way by 

i ; 
dina’ 7 
The slow-roll approximation becomes |e;| < 1 and |6;| < 
1. 

The definition of the Hubble and GB flow parameters 
renders significant simplification in the involved expres- 
sions. From Eqs. (2-4) we can express the kinetic term 
and the potential in terms of the flow parameters: 


wo? = [Qe — ô (1 + 1 — b)|H”, 
V= Le 2€4 } ô1( 5+ E1 62)|H?. 


6, = 4€H, 


di+1 = (10) 


(11) 
(12) 


We see that the potential energy dominates over the ki- 
netic energy and the GB energy. During slow roll the 
sign of w is determined by the sign of (2€, — 61). In the 
special case of 2€, = 61, the field is frozen, which corre- 
sponds to the constant Hubble parameter. We will not 
consider this special case further. 

It is known that slow roll is an attractor that is rapidly 
approached by different initial conditions [14]. Let us 
demonstrate that the slow-roll solution (5-7) is the at- 
tractor of the system (2-4) under the slow-roll condition. 
From Eggs. (3) and (4) one has 


2u (1+ 24wl gE p H4 — 4E pu H) Hy = 
—wu? + AE pou” H? 
—4t gH? (4uH + WV g +12wt 4H"), 
u (1+ 24w pE p H4 — 4E guU H) ug = 
—3uH — w (1 — 4€ gu H) Vg + 6E gH? (34? 


(13) 


—4wE gpu” H? — 2wH? + 12w gu H’), (14) 
subject to the Friedmann constraint equation 
6H? = wu? + 2V + 24€ pu H’, (15) 


where u = ¢. Suppose H(@) and u(¢) is the slow-roll 
solution to the system (13-15). Add to this a linear ho- 
mogeneous perturbation 6H(¢) and du(¢); the attrac- 
tor condition will be satisfied if it becomes small as the 
Universe expands. Inserting H(¢) = H(¢) + 6H(¢) and 
u(d) = ülo) + du(¢) into Eqs. (13-15), we find that the 
linear perturbations satisfy 


— 3H O1€4 
ôH = ae ft + = + Oet) d6H(16) 
Sug = 3H [a 4 Dense Serbs = F150 — 867 
eS u 6(2e1 — 61) 


+0(6re1,614) ou, (17) 
which have an approximately decaying solution with 
ôH x exp(—3N) and du x exp(—3N) if the Hubble and 
GB flow parameters vary slowly, and hence all linear per- 
turbations die away exponentially fast as the number of 
e-folds increases. 


III. POWER SPECTRA 


At linear order in perturbation theory, the Fourier 
modes of curvature perturbations satisfy [15] 


LA 
v” + (aw — 2) v=0, (18) 
ZR 


where a prime represents a derivative with respect to 

conformal time t = f a~'dt, and where zr and cr are 

given by 

2 a? (wd? + 6AEH?) 

R o (1— 4A)?H? ’ 

x 8A¿HH +24? H? (Ë — ÈH) 
wo? + 6AEHS 


z 


(19) 


c= i (20) 


with A = 4€H/(1 — 4€H). One can express zp and ch 
in terms of the Hubble and GB flow parameters, 


F 
2 aD 
ZR =a (@— 4A)?’ 
2e, + $61(1 — 5ey — 62) 
F ri 


(21) 


ce =1- A? (22) 
where A = 61/(1—61) and F = 21-61 (1+€1—62)+ 3.6). 
The effective mass term in the scalar mode equation (18) 
reads 
zh SaS 3 F 3 A 
ÎR = a@H?|2— a e 
mo l atar omie A 
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with 


= €1€2(2— 61) — ô1ô2(1 + €1 — 62 — ĝ3) 
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The Fourier modes of tensor perturbations satisfy [15] 


u” + (Gu = =) u=0, (24) 
ZT 
where 
22 = a*(1 — 4€H), (25) 
2 _,_ 4- EH) 
cp=1 iz XH ; (26) 


Note that the coupling £ appears not only in the k? term 
responsible for subhorizon oscillations but also in the ef- 
fective mass term z4 /zr. This differs from k-inflation in 
which the equations of motion and evolution of the ten- 
sor perturbations are not affected by nonminimal kinetic 
terms. In terms of the Hubble and GB flow parameters 
z2 and c} can be written as 


z% = a*(1 — 61), (27) 
@.=1+A( — 4 — ds). (28) 


The effective mass term in the tensor mode equation (24) 
reads 


Zp oe 3 1 
— = a‘H*|2 Ey ghee 5 Adal €1 + d2 + 63) 


ZT 


-70a . (29) 


If both €, and 6; are constants, which corresponds to 
the power-law inflation with an exponential potential and 
an exponential GB coupling [11], then Eqs. (23) and (29) 
become 

te 2p _ 1 2-4 (30) 
ZR ZT T2 (1 — €61)? 
The spectral indices of scalar and tensor perturbations 
read exactly 


2€1 
1—€’ 


nr—l=nr=-— (31) 


which are consistent with the results in Ref. [11]. We 
note that only €ı appears in the spectral indices whether 
for potential-dominated or GB-dominated inflation. 

In general the Hubble and GB flow parameters are 
functions of cosmic time. We shall assume that time 
derivatives of the flow parameters can be neglected dur- 
ing slow-roll inflation, which will allow us to obtain the 
leading contribution to the slow-roll approximation. Un- 
der this assumption one has t~! ~ —aH(1 — e1) and 
T?’zh|zr = vz — 1/4 can be approximated to be con- 
stant. Then the general solution to Eq. (18) is a linear 
combination of Hankel functions 


R 


= VAT (20m) /4 [er HH (crklrl) + cH) (erk|r|) |(32) 


We choose cy = 1 and co = 0, so that the usual 
Minkowski vacuum state is recovered in the asymptotic 
past (crk|r| + oo). The power spectrum of curvature 
perturbations Pr = k?|v/zp|?/2n? on the large scales 
(crk « aH) is 


a H (1-4/2)? 
PR = — > —__ ——— 
|F| 472? \ aHjr| 


-o 


oR 3e R H? T2 
ai 5 E ve) (33) 
|F| 4r? T? (3/2) crk= aH 
with spectral index 
NR — 1=3- 2UR. (34) 


As in the case of scalar perturbations, the power spec- 
trum of tensor perturbations Pr = 2k3|2u/zp|?/27? 
given by 


eae ae ae 
Pr = — 
1— 6; 4n? \aH |r| 
7-3 H? T? (vr) 
T Tr? T2(3/2) 


x 2Te (35) 


crk-a H 
with spectral index 
nT = 3— 20T, (36) 


where we have defined v2 = 7?z/,/zr + 1/4. All back- 
ground quantities above are evaluated at the moment 
such that crk = aH. This is not exactly the same time 
as the horizon-crossing time in Eq. (33) for scalar modes, 
but to lowest order in the slow-roll parameters this differ- 
ence is unimportant. We can use the slow-roll approxi- 
mation to estimate the amount of e-folds between horizon 
crossing of the scalar mode and the tensor mode with a 
reference scale k, AN ~ In(er/cr) ~ 61/2. An impor- 
tant observational quantity is the tensor-to-scalar ratio 
which is defined as 


(oo 


Pr 3 hT? — 
~? +2vTr—2vR F| R 


To first order in the slow-roll approximation, we have 


0? (4e1 + 61) 


2~1- 
CR 2(2e1 — 51) (38) 
ZR 2772 3(2€1€2 — 6169) 
FR — a? H = pack Sea a 
ZR . E = H 2(2e1 — 61) 
+O(ere2,518)| 5 (39) 
G~1t+h, (40) 
1 
TT = a? H? [2 — e1 + O(6162)]. (41) 


ZT 
The spectral indices of scalar and tensor perturbations 
read 
2€1€2 = ô 62 
—12—2e, — ——___ 42 
NR €1 261 — ôi ’ ( ) 
nr = —26, (43) 


which show that the spectral index of scalar perturbation 
contains not only the Hubble flow parameters but also the 
GB flow parameters. Even for a solution very close to de 
Sitter inflation (i.e., «; ~ 0), the GB term can lead to a 
red (d2 > 0) or blue (62 < 0) power spectrum of scalar 
perturbation. If |d1| < «1, the spectral indices are the 
same as for a potential-driven slow-roll inflation. 
The tensor-to-scalar ratio (37) is approximately 


r ~ 8|2e, — 61| A —8nr, (44) 


which is the modified consistency relation. The degen- 
eracy of standard consistency relation is broken in the 
slow-roll inflation with the Gauss-Bonnet correction. For 
this reason, the future experimental checking of this re- 
lation is usually regarded as an important test of the 
simplest forms of inflation. 

The Hubble and GB flow parameters can be expressed 
in terms of the potential and the GB coupling 


a x A (45) 
~ oc | Yee _ Ye , Qe 

a= -o (4 te 2e), w 

by = E sQV, (47) 
~ _p[£ss _ Ve , Qe 

by ~ o($ +e y le), (48) 


where Q = w(V¢/V + 4E V/3). 

The key result of our paper is the general slow-roll ex- 
pression for GB inflation, Eqs. (42), (43) and (44), which 
is new and follows from a nontrivial calculation. 


IV. AN EXAMPLE MODEL 


Let us consider a specific inflation model 


V(b) = Vog”, €(6) = E09". (49) 


This potential has been widely studied. The specific 
choice of GB coupling allows us to find an analytic re- 
lation between the spectral index of curvature perturba- 
tions and the tensor-to-scalar ratio. If a = 4Voéo/3 = 1, 
all flow parameters vanish. The motion of the inflaton is 
frozen because the force due to the slope of the potential 
is exactly balanced by one, the slope of the GB coupling. 
In this case, exact de Sitter inflation can be realized for 
the monomial potential and the inverse monomial GB 
coupling. If a < 1, choosing w = 1 is required for a 
positive €,. In this case the contribution of the positive 
GB term increases the Hubble expansion rate during in- 
flation, which makes the evolution of the inflaton slower 
than in the case of standard slow-roll inflation, while the 
contribution of the negative GB term decreases the Hub- 
ble expansion rate. If a > 1, we choose w = —1 to guar- 
antee €; > 0. The potential force drives the inflaton to 
climb up the potential while the GB force drives the field 
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FIG. 1: Two-dimensional joint marginalized constraint (68% 
and 95% confidence level) on the scalar spectral index nr 
and the tensor-to-scalar ratio r derived from the data combi- 
nation of WMAP7+BAO+4Hp by imposing the standard con- 
sistency relation. The symbols show the predictions from the 
$'-potential (solid line) and ¢?-potential (dashed line) models 
with the number of e-folds equal to 50 (small) and 60 (large). 


to roll down. Since the GB force dominates over the po- 

tential force, slow-roll inflation can be realized. In what 

follows we restrict our discussion to the case of a < 1. 
The flow parameters are 


a x sm —a)o-?, (50) 
e ~ 2n(l—a)d~?, (51) 
6, ~ nall- a)’, (52) 
6g ~ 2n(1-— a)’. (53) 
From Eqs. (42) and (44) one gets 
nr—-1 = —n(n+2)(1—a)d, (54) 
r = 8n7?(1—a)*¢-?. (55) 


Inflation ends at €1(¢ena) = 1, which gives the value of 
the field at the end of inflation 


Baa = En?(1 — a). (56) 


Then from (8) we find the value of the field N e-folds 
before the end of inflation 


¢? = 2n(1-a)(N + =): (57) 


The spectral index nz and the tensor-to-scalar ratio r 
can be written in terms of the function of N: 


2(n + 2) 
slas Lier A 
me AN +n’ ae) 
_ 16n(1 — a) 
r= "AN oe (59) 


Note that the spectral index is independent of Vo and £p, 
but the tensor-to-scalar ratio depends on œ = 4Voéo/3. 
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FIG. 2: Tensor-to-scalar ratio r versus the spectral index nr 
for the inflation model (49) with n = 2 (top panel) and n = 4 
(bottom panel). The contours show the 68% and 95% con- 
fidence level derived from WMAP7+BAO+Ho without the 
consistency relation. 


The GB correction leads to a reduction of the tensor-to- 
scalar ratio if éo > 0 while an enhancement if éo < 0, 
which is still valid in the power-law inflation model with 
the exponential potential and GB coupling [11]. 

Figure 1 shows the two-dimensional joint marginalized 
constraint (68% and 95% confidence level) on nr and 
r from the 7-year WMAP+BAO+Ho by imposing the 
standard consistency relation [16]. The symbols show 
the predictions from the ¢4-potential (solid line) and ¢?- 
potential (dashed line) models with the number of e-folds 
equal to 50 (small) and 60 (large). We can see that the 
predicted points with N = 50,60 for the quartic poten- 
tial are far away from the 95% region. The quadratic 
potential is consistent with the data. 

However, the consistency relation nr = —r/8 is bro- 
ken in the slow-roll inflation with the GB correction. 
Therefore, in our analysis, ny is varied independent of 
the tensor-to-scalar ratio. For the tensor perturbations 
we assume a power-law power spectrum, with a uniform 
prior on np as —0.5 < nr < 0. In Fig. 2 we show 
the lo and 20 contours derived from the data combi- 
nation of WMAP7+BAO+Hp by using the CosmoMC 
package [17]. Compared to the contours of Fig. 1 we find 


that the joint constraint on nr and r becomes a little 
tighter. The WMAP7+BAO+4p data do not constrain 
nr. Basically all values allowed by the prior are also 
allowed by the potential and the coupling. 

In Fig. 2 we plot the values of nr and r in the models 
with n = 2 (top panel) and n = 4 (bottom panel) for 
different values of N and a. We can see that the model 
parameter a can shift the predicted r vertically for a fixed 
number of e-folds. For n = 2, the model with a positive 
a is more favored observationally. For n = 4, the model 
with a > 0.7 is consistent with the data within the 95% 
confidence level, in which the prediction for the tensor- 
to-scalar ratio is smaller than the a = 0 case while the 
prediction for nz is the same as the a = 0 case. Other 
ways to avoid the exclusion of the ¢* potential have been 
studied in Ref. [18]. 


V. CONCLUSIONS AND DISCUSSIONS 


In this paper we have studied slow-roll inflation with 
a nonminimally coupled Gauss-Bonnet term. We have 
defined a combined hierarchy (¢;,6;) of Hubble and GB 
flow functions such that |e;| < 1 and |é;| < 1 is the 
analogue of the standard slow-roll approximation. It has 
been demonstrated that slow-roll solution is the attrac- 
tor solution under the slow-roll condition. We have an- 
alytically derived the power spectra of scalar and tensor 
perturbations. In general the spectral index of scalar per- 
turbations depends on the Hubble flow parameters and 
the GB flow parameters. However, the spectral index of 
tensor perturbations is independent of the GB flow pa- 
rameters to first order in the slow-roll approximation. In 
this scenario the standard consistency relation does not 
hold because of the GB correction. 

We apply our general formalism to large-field inflation 
with a monomial potential and the GB coupling (49). 
We focus on the case of w = 1 and a < 1 since the 
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field theory of phantom-type fields encounters the prob- 
lem of stability. In this case, the GB term with the pos- 
itive (or negative) coupling slows down (or speeds up) 
the evolution of the inflaton during inflation, which de- 
creases (or increases) the energy scale of the potential 
to be in agreement with the amplitude of scalar pertur- 
bations. However the amplitude of tensor perturbations 
only depends on the energy scale of the potential at the 
horizon-crossing time. Therefore, the tensor-to-scalar ra- 
tio is suppressed for œa > 0 while it is enhanced for a < 0. 

As shown in Fig. 2, the model parameter & can shift 
the predicted r vertically for a fixed number of e-folds 
in the nr-r plane. For n = 2, the quadratic potential 
can be made a better fit to the data by the positive GB 
coupling. For n = 4, it is known that the model with 
a = 0 is excluded by the WMAP7+BAO+AHp analysis. 
However, in our scenario of inflation a > 0.7 is within 
the 20 contour for N > 50, and it is consistent with the 
data within the 95% confidence level. 

The results of this work are generic as soon as nonmin- 
imal couplings are considered. While it is always possible 
by means of a conformal transformation to work in the 
Einstein frame and to avoid the presence of a ¢?R term 
in the Lagrangian, the coupling of the scalar field to the 
GB term cannot be argued away by the same conformal 
transformation. While we studied perturbation spectra 
in the Einstein frame, similar properties hold in the Jor- 
dan frame. 
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